Does interferometry probe thermalization? 



On 
O 
O 
(N 

>v 



Clement Gombeaud, 1 Tuomas Lappi, 1,2 and Jean- Yves Ollitrault 1 

'Institut de Physique Theonque, CEA/DSM/IPhT, CNRS/MPPU/URA2306 

CEA Saclay, F-91191 Gif-sur-Yvette Cedex. 
2 Department of Physics P.O. Box 35, 40014 University of Jyvaskyla, Finland 

(Dated: May 26, 2009) 

We carry out a systematic study of interferometry radii in ultrarelativistic heavy-ion collisions 
within a two-dimensional transport model. We compute the transverse radii R and R s as a function 
of pt for various values of the Knudsen number, which measures the degree of thermalization in the 
system. They converge to the hydrodynamical limit much more slowly (by a factor ~ 3) than elliptic 
flow. This solves most of the HBT puzzle for central collisions: Ro/Rs is in the range 1.1 — 1.2 for 
realistic values of the Knudsen number, much closer to experimental data (~ 1) than the value 1.5 
from hydrodynamical calculations. The pt dependence of Ro and R s , which is usually said to reflect 
collective flow, also has a very limited sensitivity to the degree of thermalization. We then study the 
azimuthal oscillations of R , R a , and R os for non central collisions. Their amplitudes depend little 
on the Knudsen number, and reflect the eccentricity of the overlap area between the two nuclei. 

PACS numbers: 25.75.Gz 
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I. INTRODUCTION 

Correlations of identical particles produced in ultrarel- 
ativistic heavy-ion collisions have the unique capability 
to access directly the size of the fireball More pre- 
cisely, they measure the separation distribution of par- 
ticles with a given momentum p (regions of homogene- 
ity (H) after the last interaction. These data, referred 
to as HBT Q, thus impose severe constraints on model 
calculations. In particular blast-wave 0] and hydrody- 
namical models [5j, la, [3,18, 9], which have been rather suc- 
cessful in reproducing transverse momentum spectra and 
elliptic flows of identified particles up to p t — 2 GeV/c, 
fail in reproducing HBT radii. More specifically, they 
generally overpredict the longitudinal size Ri , as well as 
the ratio R Q / R s , where R D and R s are the transverse radii 
parallel and orthogonal to the transverse momentum, re- 
spectively. On the other hand, they correctly predict the 
decrease of radii with p tl which is often claimed to be 
a signature of collective flow. Viscous hydrodynamics 
gives smaller values of R /R s than ideal hydrodynam- 
ics (THE)]. Transport models also yield a smaller value, 
much closer to data, typically around 1.2 [12, El, Hi- 

In this paper, we investigate systematically the sensi- 
tivity of HBT radii to the degree of thermalization in the 
system. We explain the difference between predictions 
from hydrodynamics (where local thermalization is as- 
sumed) and transport models (where the system is gen- 
erally not locally equilibrated). We consider a simple 
model, where the system consists of massless particles 
undergoing 2 — > 2 elastic collisions [151 ] . The mean free 
path of a particle between two collisions can be chosen 
arbitrarily by varying the cross section. The limit of zero 
mean free path is the "hydrodynamic limit" : the sys- 
tem is locally thermalized and its expansion follows the 
laws of ideal hydrodynamics. The limit of infinite mean 
free path corresponds to free-streaming particles: in this 
case, HBT radii reflect the initial distribution of parti- 



cles. For finite values of the mean free path, the system 
is partially thermalized. In this paper, we study quanti- 
tatively how HBT observables vary between these two ex- 
tremes. A further simplification is that we consider only 
a two-dimensional system living in the transverse plane: 
the present study is therefore limited to the transverse 
radii R Q and R s (and the cross term R os for non-central 
collisions), and longitudinal expansion is not taken into 
account. 

We do not mean to provide a realistic model of heavy- 
ion collisions. The fact that we consider only 2^2 
collisions implies, for sake of consistency, that the sys- 
tem is dilute (higher-order processes such as 3 — > 3, are 
negligible), hence the equation of state is that of a per- 
fect gas. This is the price to pay for a control handle on 
thermalization. In the real world, the equation of state 
of QCD is not that of a perfect gas: it has a very sharp 
structure around T ~ 170 MeV [la], which is expected to 
influence observables, including HBT radii. By compar- 
ing with experimental data from RHIC, we expect to fail 
whenever the equation of state is important. This will 
allow us to disentangle effects which can be attributed 
to the equation of state, from those which are due to 
thermalization and flow. 

This article is organized as follows. In Sec. HH we 
present our model and explain how HBT radii are ob- 
tained. Sec. IIIII discusses the pt dependence of R and 
R s and the value of R /R s for central collisions. Results 
from transport theory and ideal hydrodynamics are com- 
pared. Sec. |IV] discusses the azimuthal oscillations of 
R s and R os for noncentral collisions. Our conclusions are 
summarized in Sec. [V] 



II. MODEL 

Nuclei colliding at RHIC are thin pancakes due to the 
strong Lorentz contraction along the collision axis. This 
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large separation between the longitudinal and transverse 
scales implies that longitudinal and transverse dynamics 
are to a large extent decoupled. In this paper, we concen- 
trate on the transverse expansion, which we model using 
a 2-dimensional relativistic Boltzmann equation [15j . We 
first describe the initial conditions of the evolution. We 
briefly recall how the Boltzmann equation is solved. We 
then define the Knudsen number, which measures how 
close the system is to local thermal equilibrium. We fi- 
nally define HBT radii. 

A. Initial conditions 

The nucleus-nucleus collision creates particles. We as- 
sume for simplicity that the spatial distribution of these 
particles is initially a gaussian in the transverse plane: 

n(x,y) = - e ^ , (1) 

A-KOxOy 

where N is the total number of particles, and a x and a y 
are the rms widths of the distributions in the x and y 
directions. The x axis denotes the direction of impact 
parameter, or reaction plane. 

As for the initial momentum distribution, two differ- 
ent scenarios have been implemented and compared. The 
first scenario is the same as in (l5j . In order to compare 
transport theory and hydrodynamics, we take the same 
initial conditions: The momentum distribution is locally 
thermal, and the temperature is related to the density 
according to the equation of state of a 2-dimensional 
massless, ideal gas: T oc n 1 / 2 . Since our calculation is 
purely classical, we assume Maxwell-Boltzmann statistics 
for sake of consistency: 

where T{x, y) is the local temperature, given by: 

T(x ' y)=ToeXP (~^I"^f)- (3) 

The second set of initial conditions are taken from 
the Color glass condensate (CGC) calculations (l7l [T^|. 
where the initial gluon spectrum is calculated by solving 
the classical Yang-Mills (CYM) equations with the initial 
conditions given by the MV model [l9| . The result of the 
numerical computation can be parameterized 20] as 

dN f ai [ e &_i]-i (p t /A s )<1.5 
cPp t (Px \ a 2 log(4^ t /A s )(p t /A s )- 4 {p t /A s ) > 1.5 

(4) 

with on = 0.137, a 2 = 0.0087 and b = 0.465. The color 
charge density parameter A s (proportional to the satu- 
ration scale Q s [2l|) plays the role of the temperature 
as the only transverse momentum scale in the system. 
The parameterization Q was fit to a calculation for a 



nucleus of an infinite size on the transverse plane, but we 
generalize it by letting A s (x,y) have the same Gaussian 
dependence on the transverse coordinate as the temper- 
ature in Eq. ([3]), with an absolute value adjusted to give 
the same value for {pt}- 

Our 2-dimensional kinetic theory approach does not 
contain longitudinal expansion, and therefore cannot ad- 
dress questions related to isotropization of the particle 
distribution. The CGC initial conditions naturally lead 
to a very anisotropic initial condition where, after r ~ 
1/Qs) (Pz) •C (j>t), whereas conventional 3 dimensional 
hydrodynamics assumes isotropy in the local rest frame. 
In the 2-dimensional approach, p z = for all particles, 
and the energy per particle {pt} is constant throughout 
the evolution; we adjust it to (pt) = 420 MeV, corre- 
sponding roughly to the value for pions at the top RHIC 
energy J22||. This fixes the value of To for the thermal 
initial conditions Q {(pt) = f^o) an d the value of A s 
in Eq. (01 for the CGC initial conditions. In practice, 
this implies an unrealistically small value of the satura- 
tion scale Q s . Conventional estimates of Q s are larger, 
but significant longitudinal cooling is required in order 
to match observed pt spectra. In our calculation, the 
conservation of {pt) makes our initial {pt} smaller than 
most estimates. However, our emphasis in this paper is 
on the influence of thermalization on the transverse HBT 
radii and the transverse momentum dependence. For this 
purpose the two initial conditions that we use represent 
the opposite ends of the range of physically reasonable pt 
spectra: from fully thermalized to one with a perturba- 
tive power law behavior at large pt that one would expect 
to observe in the absence of any final state interactions. 

B. Expansion: Knudsen number 

The results presented in this paper use the algorithm 
described in [l5| to solve the two-dimensional relativistic 
Boltzmann equation. The Monte-Carlo algorithm follows 
the trajectory of every particle throughout the expansion 
of the system, until they cease to interact. Particles in- 
teract through 2^2 elastic collisions. The cross section 
is assumed isotropic in the center-of-mass frame for sim- 
plicity. 

The only remaining parameters in the simulation are 
the total number of particles, N, and the elastic cross 
section, a. (Since we are working in two dimensions, 
a has the dimension of a length.) The initial average 
particle density per unit surface is: 



It is worth emphasizing that N is usually much larger 
in the Monte-Carlo simulation than in an actual heavy- 
ion collision ("parton subdivision" technique). This can 
be understood in the following way. The physical length 
scale that should be independent of N is the mean free 
path A. Its precise value depends on the initial velocity 
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and position of the particle, but the order of magnitude 
is generally A = l/an. A must be compared to another 
length scale, the average interparticle distance d — nT 1 ! 2 . 
Let us define the dilution parameter D as 
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(6) 



Our description of the system in terms of elastic 2^2 
collisions is consistent only in the limit when D is small 
and the contribution of many-body collisions is sup- 
pressed. This is a requirement of the Boltzmann equa- 
tion [l5| ■ To achieve this one must take the limit of large 
N and small a keeping aN fixed. This ensures that our 
results are extrapolations to the limit N — > oo at fixed A. 
Because in this limit a approaches zero the interactions 
between the particles become truly pointlikc, and prob- 
lems with causality and Lorentz-invariance are avoided. 
For this reason, all the results presented in this paper are 
obtained by doing two simulations with the same value 
of A and different values of D\ the results are then ex- 
trapolated linearly to D = 0. 

The standard dimensionless parameter to characterize 
the degree of thermalization is the Knudsen number K, 
defined as the ratio of the mean free path to the charac- 
teristic size of the system R. We define R as in [23| : 
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The Knudsen number K is then defined as 
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The inverse of the Knudsen number is proportional to 
the average number of collisions per particle, n co \\. The 
product n co \\K remains very close to 1.6, for all values 
of K [15(. Hydrodynamics is the limit K — > 0, while 
K — > +oo correspond to free streaming particles. A fit 
to the centrality dependence of elliptic flow [2J| suggests 
that K ~ 0.3 for central Au-Au collisions at RHIC. 

We choose to keep the scattering cross section a con- 
stant as a function of time for sake of simplicity, as for 
instance in the AMPT transport model [25|. Other trans- 
port calculations have been carried out |26| where the 
viscosity to entropy ratio r]/s is kept constant, so that a 
depends on temperature or time, typically like i 2 / 3 . As 
we recall below, HBT radii give a measure of the sys- 
tem when the last scattering occurs, that is, much later 
than other observables such as elliptic flow. Therefore, 
results might differ significantly with a time-dependent 
cross section. 



C. HBT radii 

For a particle with momentum p t , we denote by 
(t, x, y) the space-time point where the last collision oc- 
curs. The "out" and "side" coordinates are then defined 



as the projections parallel and orthogonal to the particle 
momentum: 

x a = x • v — vt = x cos 4> + y sin 4> — vt 

x s = x x v = x sin (f> — y cos 0, (9) 

where v = pt/pt is the particle velocity (v = 1), and <p 
its azimuthal angle: p t = (p t cos (f>,pt sin</>). Both x Q and 
x s are invariant under a translation along the trajectory 
after the last scattering: (t, x] — > (t + r, x + vt) . In 
particular, they are invariant through a scattering at zero 
angle. HBT radii are obtained by averaging over many 
particles with the same momentum: 
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(10) 



Radii defined in this way coincide with those obtained 
from the curvature of the correlation function at zero 
relative momentum, in the absence of final-state inter- 
actions [l2 |. Experimentally, radii are usually obtained 
from gaussian fits to the correlation function. This proce- 
dure gives different radii if the source is not gaussian [l[ . 
In our case, we have checked explicitly that sources are 
close to gaussian, but we have not investigated system- 
atically effects of non gaussianities. Strictly speaking, 
averages in Eq. (|10p are for a given momentum. In prac- 
tice, our results are obtained by taking bins of width 
10 MeV/c in p t) and averaging over the particles in the 
bin. We have checked that results do not vary signifi- 
cantly with a smaller bin size. 

The radii defined by Eq. (fT0|) are generally functions 
of pt and </>. In Sec. IIII| we study central collisions, with 
'x = a y Symmetry of the system with respect to the di- 
rection of p t then implies R os = 0. Rotational symmetry 
implies that R and R s are independent of <p. The more 
general case when a x and a y differ is studied in Sec. IIVI 



III. CENTRAL COLLISIONS 

In this section, we discuss how the pt dependence of R Q 
evolves with the Knudsen number for a central collision. 
We then discuss the ratio R /R s - Finally, we compare 
with experimental data. In order to mimic a central Au- 
Au collision at RHIC, we use the initial density profile 
Eq. (fTJ) with o x = a v = 3 fm. This value corresponds to 
the rms width of the initial density profile in an optical 
Glauber calculation (27j . 

The decrease of HBT radii with the transverse momen- 
tum p tl typically like p^ 1 ^ 2 ' , is often [H, [H, [29[ presented 
as a signature of collective flow. Collective flow is as- 
sociated with the hydrodynamic limit, i.e., the limit of 
small K. Fig. [T] displays R D versus pt for thermal ini- 
tial conditions, Eq. ([2]), and several values of the Knud- 
sen number K . Generally, R increases as K decreases. 
However, this is a small effect. The decrease of R with 
p t is more pronounced in the hydrodynamic limit (small 
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FIG. 1: HBT radius R versus transverse momentum pt of 
particles in the transport calculation. The curves are labeled 
by the value of the Knudsen number K. 
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FIG. 2: HBT radius R versus transverse momentum p t of 
particles in ideal hydrodynamics, at a given time t. The curves 
are labeled by the value of t. 



K) but is also seen for free streaming particles (large 
K). For large K, HBT radii reflect the initial momen- 
tum distribution: both with thermal initial conditions 
conditions, Eq. ([2]), and with CGC initial conditions, 
Eq. (HJ), particles with higher p t are more likely to be 
produced in dense regions, i.e., near the center of the 
fireball x = y = 0. For p t > (p t ), Eqs. © and © 
yield R {pt) — o- x y / 1.5{pt) /pt for the initial distribution, 
while Eq. ((H) gives R (pt) — o x /^/2. In practice, after 
collisions have occurred, both sets of initial conditions 
yield similar radii, as we shall see explicitly later. 

Decreasing the Knudsen number K amounts to in- 
creasing the number of collisions, hence the "freeze-out" 
time when the last collision occurs. To show this ex- 
plicitly, Fig. [5] displays R (pt) in ideal hydrodynamics, 
assuming sudden freeze-out at time t. The equations of 
hydrodynamics are solved using a first-order Godunov 



FIG. 3: R and R s , averaged over the interval 0.25 < p t < 
0.75 GeV/c, versus 1/K, which scales like the number of 
collisions per particle. The lines are 3-parameter fits with 
Eq. (HIJ. The dotted curve shows, for sake of illustration, the 
variation of elliptic flow in a non-central collision, scaled by 
the hydrodynamical limit (from [15J). 



scheme [30j. For sake of consistency with the transport 
calculation, the equation of state of the fluid is that of 
a two-dimensional ideal gas, and there is no longitudinal 
expansion [15j . Hydrodynamics at t — 0.2 fm/c gives 
the same radii as the transport calculation for large K in 
Fig. [U because we have chosen the same initial conditions 
for both calculations. As time evolves, R increases; the 
increase is more pronounced and occurs later at low p t . 
At a given p t , the value of R converges as t increases. 
This is by no means a trivial result: the location of the 
last interaction, (x Q ) in Eq. (|10[) . increases linearly with 
t. Only the dispersion R Q of this location converges. Hy- 
drodynamics at large t is almost identical to transport 
at small K (the relative difference is less than 5%). This 
is also a non-trivial result, although it is implicit in all 
hydrodynamical studies of HBT observables [5j: HBT 
observables are defined at the last scattering, when the 
system is no longer in local equilibrium, and hydrody- 
namics is not valid. 

Hydrodynamical calculations usually yield a value of 
R /R s which is much too large, of the order of 1.5, while 
RHIC data are compatible with 1. While R Q increases 
with time in hydrodynamics, R s decreases. Initially, 
Ro{vt) = Rs{pt) by symmetry. In the transport calcu- 
lation, the same behavior is observed as the number of 
collisions per particle 1/K increases, as shown in Fig. [3] 
However, this increase is quite slow. The same curve 
shows, for sake of illustration, the increase of the elliptic 
flow i>2 with 1/K for a noncentral collision 15j . Elliptic 
flow converges to the "hydrodynamic limit" much faster 
than HBT radii. In order to put this statement on a 
quantitative basis, we fit our numerical results for R a (K) 
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FIG. 4: R /R s versus p t for AT = 0.3. Dashed line: thermal 
initial conditions, Eq. ((2|; full line: CGC initial conditions, 
Eq. 0. Error bars on the transport calculations are statisti- 
cal only. Stars: data from STAR (3l|. 



1.5. It has been recently argued that ideal hydrodynam- 
ics with an early freeze-out [13] also explains the HBT 
puzzle. Our results also show that increasing the Knud- 
sen number amounts to decreasing the freeze-out time in 
ideal hydrodynamics. Viscous corrections to ideal hydro- 
dynamics, which incorporate deviations to local equilib- 
rium to first order in the Knudsen number K, also lead 
to a reduced R /R s (together with a reduced longitudi- 
nal radius Rl, also closer to data). However, the value 
of rj/s required to match [To| the data is much larger 
than that inferred from the study of elliptic flow [33| , so 
that viscous corrections explain only a small part of the 
HBT puzzle [ll|. However, viscous hydrodynamics itself 
breaks down at freeze-out, and may not be a reliable tool 
for estimating HBT radii. Our results suggest that devi- 
ations from equilibrium have larger effects on HBT radii 
than inferred from viscous hydrodynamics. We find that 
partial thermalization, which has been shown to explain 
the centrality dependence of v 2 , also solves most of the 
HBT puzzle for R a /R s . 



with the following formula [If 



R (K) = Rl' 



jjhydro _ R f, 

1 + K/K 



(11) 



The fit parameters are the free-streaming (K — > oo) limit 
R[ s -, the hydrodynamic (K -> 0) limit R^ dl °, and K , 
the value of the Knudsen number for which R (K) is half- 
way between free-streaming and hydro. A similar formula 
can be used for R S (K). It fits our numerical results per- 
fectly (see Fig. EJ| . The value of K is 0.167 ± 0.007 for 
R Q and 0.215 ± 0.006 for R s , while it is 0.7 for v 2 HI: 
convergence toward the hydrodynamic limit requires 3-4 
times more collisions for HBT radii than for elliptic flow. 
The other fit parameters are R% ydl ° = 3.071 ± 0.006 fm 
and R^ dl '° = 2.088 ± 0.006 fm: we recover the HBT 
puzzle R I R s ~ 1.5 in the hydrodynamical limit. 

We now discuss the value of R /R s at RHIC. The cen- 
trality dependence of vi suggests that K ~ 0.3 in cen- 
tral Au-Au collisions (24j. For this value, v-i is already 
70% of the hydrodynamic limit. On the other hand, 
R /R s ~ 1.16, which is significantly below the hydro- 
dynamic limit of 1.5. A similar value (R n / R s ~ 1.2) 
was found with the AMPT transport code [12|]. A more 
detailed comparison is shown in Fig. [H which displays 
R /R s versus pt for the two sets of initial conditions, 
Eqs. and (j4|). K has been fixed to the value which is 
favored by v-i data [HI], i- e -> K = 0.3 for central collisions. 
For both sets of initial conditions, R a /R s is essentially 
independent of pt, while data show a slight decrease. In 
this respect, our results differ from the covariant MPC 
model (which is in principle equivalent to ours, with the 
longitudinal expansion taken into account), where R a / R s 
is found smaller than 1 at large p t [IH • 

Our value of R /R s for central collisions is in much 
better agreement with experimental data than models 
based on ideal hydrodynamics, which give a value around 
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FIG. 5: HBT radius R a versus transverse momentum pt. 
Lines: our results for K = 0.3 and different initial conditions: 
thermal, Eq. and CGC, Eq. HJ. Symbols: experimental 
data from STAR [H, PHOBOS and PHENIX JH- 



While our transport calculation gives a plausible ex- 
planation for the small R /R s , it completely misses the 
absolute value of HBT radii. This is shown in Fig. [5l 
which displays a comparison between R (p t ) fr om our 
transport calculation with data from STAR [3l|, PHO- 
BOS [H and PHENIX Experimental values are 
much larger. This is due to the equation of state @,[36|], 
which is that of an ideal gas in our calculation. Our HBT 
volume R R S is essentially independent of the number of 
collisions 1/K. It has been argued that this is a general 
result for an ideal gas, due to entropy conservation [37l ]. 
The equation of state of QCD, on the other hand, has a 
sharp structure around T c ~ 170 MeV: as the tempera- 
ture decreases, the entropy density drops by an order of 
magnitude in a narrow interval around T c . In a heavy- 
ion collision, the volume increases by a large factor with 
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essentially no change in the temperature. This explains 
why HBT volumes increase at the transition. A com- 
plete study must also take into account the longitudinal 
expansion, and its effect on the longitudinal radius R^. 



IV. AZIMUTHALLY SENSITIVE HBT 

For a non-central collision, the interaction region is el- 
liptic, and HBT radii depend on 0. Azimuthally-sensitive 
interferometry has been investigated theoretically within 
hydrodynamical models [H, [3!J H3, EH, E3] and trans- 
port models (43[. We first briefly recall why and how 
radii depend on 0. We then study how the various radii 
depend on the Knudsen number. Finally, we introduce 
dimensionless ratios of oscillation amplitudes, which do 
not seem to have not been studied previously, and we 
compare our results with experimental data [3l| . 




FIG. 6: Illustration of dependent HBT radii. 

Fig. [S] illustrates the dependence of transverse radii. 
The initial distribution of matter is elongated along the 
y axis, so that R Q has a maximum at = tt/2, while R s 
has a maximum at — 0. Finally, R os differs from zero 
when the principal axes of the region of homogeneity are 
tilted relative to the direction of momentum. 

Before we present our results, let us briefly explain how 
the dependence of the radii is evaluated in the Monte- 
Carlo solution of the Boltzmann equation. The radii (flU)) 
involve average values, such ), which depend on 0. 

Such averages can be computed by binning in 0, and 
computing the average in each bin: 



N 



(12) 



where N is the number of particles in the bin, which 
depends on due to elliptic flow. Since the dependence 
is smooth, more accurate results are obtained by fitting 
both the numerator and the denominator of Eq. (|12p by 
Fourier series, using the symmetries — » — and ~ * 
+ 7T to restrain the number of terms (44j. Since the 
Fourier expansion converges rapidly, we keep terms only 
up to order cos 40 and sin 40. 

Fig. [7] displays the dependence of R 2 , R 2 , and R os . 
The variation of R\ and R 2 S is clearly dominated by a 
cos 20 term, while the variation of R os goes like sin 20. 
The mean value of R 2 is slightly higher than the mean 
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FIG. 7: Azimuthal dependence of HBT Radii. Thermal initial 
conditions, with o x — 1.95 fm, a y = 2.6 fm, corresponding 
roughly to a Au-Au collision at RHIC with impact parameter 
b — 7 fm. The p t interval is the same as in Fig. [3] The 
Knudsen number is K = 0.4. 



value of R 2 , which is not surprising since final-state inter- 
actions increase R and decrease R s . At = 0, however, 
R < R s , reflecting the initial eccentricity of the system. 
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FIG. 8: In-plane (0 = 0) and out-of-plane (0 = n/2 radii 
versus K~ x for thermal initial conditions. Initial conditions 
and pt interval as in Fig. [7] Lines are 3-parameter fits using 
Eq. dTTJ. 

Fig. [8] displays radii in the reaction plane (0 = 0) 
and out of the reaction plane (0 = 7r/2) versus K. 
R increases and R s decreases as the number of col- 
lisions 1/K increases, as already observed for central 
collisions (Fig. [3]). Upon closer scrutiny, Fig. [5] reveals 
that the slope of the curves differ. This is reflected 
by the value of the parameter fit Kg in Eq. (fTTj) . Kq 
is largest for R o (0) (K = 0.38 ± 0.01), smallest for 
R (it/2) (Kq = 0.04 ±0.05), and intermediate for R s (0) 
and R s (n/2) (K = 0.26 ± 0.03 and 0.20 ± 0.02, respec- 
tively) . Our interpretation is that thermalization is faster 
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in plane than out of plane, which is natural since collec- 
tive flow is preferentially in plane. 

We now study quantitatively how oscillation ampli- 
tudes vary with K. There are three such amplitudes, 
as illustrated in Fig. [7] 



AR 2 Q = R 2 (n/2)-R 2 o (0) 
AR 2 S = R 2 (0)-R 2 (n/2) 
AR os = R os {3n/4)-R os (ir/4). 



(13) 



In Fig. [JJ all three amplitudes are clearly comparable. 
If K ^> 1, particles escape freely after they have been 
produced. Setting t = in Eq. ^ and using the fact 
that the initial distribution is centered at x = y = and 
has y — > —y symmetry, one easily shows that all three 
amplitudes are equal to (y 2 — a; 2 ). The results are inte- 
grated over the pt range 0.25 < pt < 0.75 GeV/c, but our 
results depend weakly on p t . In particular, we do not see 
the inversion of oscillations at large pt reported in earlier 
hydrodynamical calculations [H, |4(|. This inversion was 
not observed in more recent calculations 14211. 





b 1AK data 


our results 


K = 0.32 


K = 0.51 


AR /AR a 


1.45 =b 0.61 


1.08 ±0.02 


1.05 ±0.02 


AD /A T)^ 

AR 0S /AR 


0.68 ± 0.42 


0.97 ±0.03 


0.99 ±0.03 




U.UoU ± U.UzD 


0.205 ± 0.003 


0.213 ±0.005 




STAR data 


our results 


K = 0.31 


K = 0.49 


ARi/ARi 


1.09 ±0.46 


1.14 ±0.02 


1.06 ±0.02 


AR os /ARi 


0.65 ±0.31 


0.90 ± 0.04 


0.92 ± 0.04 




0.086 ±0.017 


0.172 ±0.005 


0.174 ±0.005 



TABLE I: Comparison between results from STAR (3l| and 
our calculations. Top: centrality interval 20-30% and kt € 
[0.15, 0.25] GeV/c. Bottom: centrality interval 10-20% and 
k t G [0.35, 0.45] GeV/c. 



seen in HBT radii, for instance in R s : 
, _R 2 (0)-R 2 (n/2) 



s " R 2 (0) 

with the initial eccentricity 



^(tt/2) 
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FIG. 9: Ratios of oscillation amplitudes versus K' 1 for ther- 
mal initial conditions. _R„, R 2 S and R oa are integrated over the 
pt interval 0.5 < pt < 0.75 GeV/c. Lines are drawn to guide 
the eye. 

Oscillation amplitudes scale like the eccentricity of the 
overlap area between the two nuclei, which depends on 
centrality and is not known directly. This dependence 
can be avoided by considering ratios of oscillation am- 
plitudes. Out of 3 amplitudes, one may construct 2 ra- 
tios, AR 0S /AR 2 and AR 2 a /AR 2 . These ratios can be 
extracted directly from experimental data, and are equal 
to unity in the free-streaming limit (large K). They are 
plotted in Fig. [9]versus 1/K. Final-state interactions in- 
crease the oscillations of R relative to R s , much in the 
same way as they increase R relative to R s for central 
collisions. The opposite behavior was found in hydro 
and blast-wave [J] calculations, and we do not understand 
the origin of this discrepancy. 

For realistic values of K, both ratios deviate little from 
unity. It is also interesting to compare the eccentricity 



In the limit K — > oo, e s and e are strictly equal for both 
sets of initial conditions. The ratio e s /e is plotted in 
Fig. [HI It also remains close to unity. We conclude that 
none of the observables we can construct from oscillation 
amplitudes is an interesting probe of thermalization. For 
realistic values of K, the region of homogeneity essen- 
tially retains the shape of the initial distribution. 

Although our model calculation is too crude to repro- 
duce the magnitude of HBT radii, we expect that the 
above ratios are less model dependent; in particular, they 
all go to 1 in the absence of final-state interactions. A 
comparison with existing data is therefore instructive. 
Table H] displays comparisons between our results and ex- 
perimental data from STAR [3l[. The correspondence 
between centrality and eccentricity was taken from [45] ]. 
For each set of data, the two values of the Knudsen num- 
ber span the range inferred from the centrality depen- 
dence of V2 [13]. Note that the pt ranges differ for the 
two centrality intervals. This is the reason why our re- 
sults are also slightly different, although the values of K 
are essentially the same. Our results for AR 2 /AR 2 and 
AR os / AR 2 are compatible with experimental data, but 
the latter have large error bars. On the other hand, the 
experimental value of e s is smaller by a factor 2 than 
our value. In our calculations, e s remains very close to 
the initial eccentricity (see Fig.[S]). Experimentally, how- 
ever, the initial eccentricity seems to be washed out by 
the expansion. This is a spectacular effect, whose im- 
portance doesn't seem to have been fully appreciated so 
far. Hydrodynamical calculations have been reported [42j 
which arc in fair agreement with the measured value of e s . 
These calculations use a soft equation of state: it is likely 
that the soft equation of state of QCD is responsible for 
the reduced eccentricity seen in data. 
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V. CONCLUSIONS 

We have carried out a systematic study of how HBT 
observables evolve with the degree of thermalization in 
the system, characterized by the Knudsen number K. 
The number of collisions per particle scales like 1/K, and 
local equilibrium corresponds to the limit K — ► 0. Our 
results show that HBT observables depend very weakly 
on K: 

• A decrease of R with p t is expected from initial 
conditions; collective flow only makes this decrease 
slightly stronger. 

• The ratio R /R s increases very slowly when one 
approaches the hydrodynamical limit. For the val- 
ues of K found in Ref. [HJ|, it is lower than 1.2, 
and much lower than predicted by hydrodynamics. 
Partial thermalization solves most of the HBT puz- 
zle. 

• For non-central collisions, the variations of R 2 , R 2 
and R os with azimuth have almost equal ampli- 
tudes. The final eccentricity seen in the side radius 



R 2 S is very close to the initial eccentricity. 

Our results are in quantitative agreement with data for 
Ro/Rs, AR 0S /AR 2 and AR 2 JAR 2 S . On the other hand, 
our absolute values for R a and R s are much too small. 
Experimentally, it is also found that the final eccentric- 
ity is smaller than the initial eccentricity, almost by a 
factor 2. Both effects cannot be due to flow alone. On 
the other hand, they might be a signature of the soft- 
ness of the QCD equation of state or, equivalently, of the 
transition from a quark-gluon plasma to a hadron gas. 
When the quark-gluon plasma transforms into hadrons, 
the volume of the system increases by a large factor: the 
source swells, which results in larger radii and a smaller 
eccentricity. 
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